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We analyze the sharpness of crossing ("isosbestic") points which are observed in many quantities 
described by a function f{x,p), where a; is a variable (e.g., the frequency) and p a parameter (e.g., 
the temperature). We derive a simple criterion that needs to be fulfilled by f{x,p) to exhibit a sharp 
isosbestic point and relate it to the applicability of perturbation theory in p. Using this approach we 
explain the sharpness of crossing points in the optical conductivity a (lu, n) of the Falicov-Kimball 
model and the spectral function A (oj, U) of the Hubbard model. We also analyze isosbestic points 
in data for the Raman response xi^jT) of HgBa2Cu04+5, photoemission spectra I {uj,T) of thin 
VO2 films, and the reflectivity R{u,T) of CaCu3Ti40i2. 

PACS numbers: 



I. INTRODUCTION 

Already in 1937 it was noticed that the curves of 
the extinction coefficient q;(cli, ci) of a mixture of two Hq- 
uid solutions intersect exactly at one frequency w* when 
plotted as a function of uj for different concentrations ci 
of one of the components. At this particular frequency 
the extinction coefficient therefore no longer depends on 
the parameter ci . This point of intersection of the family 
of curves a(a;,ci) is referred to as isosbestic point |2H4]. 
Today isosbestic points of the extinction coefficient of two 
compounds play an important role in spectrophotometry 
[5. 

It is a common practice in the natural sciences to il- 
lustrate the functional dependence of a quantity / on a 
variable x and a parameter p by plotting f{x,p) vs. x for 
several values of p, say pi, . . .p„. This leads to a family 
of n curves f(x,pi),i — l,...n. The curves will usually 
intersect. If the points of intersection are confined to a 
narrow region of x values, or if they even coincide at a 
single point x* as in the case of the absorbance a(a;, Ci), 
it leads to a conspicuous feature whose origin calls for 
an explanation. Well-defined crossing points of curves 
are now often called "isosbestic points" even if they do 
not fall onto a single point |6I7| . We will use this term 
also in the present paper. We stress that the characteris- 
tic p-independence of f{x,p) at x* is these general cases 
only fulfilled in a fairly small p-interval, i.e., it is a "local" 
phenomenon. 

For example, an unexpected sharp crossing point at 
T* = 160mK was recognized in the family of curves of 
the specific heat cy (T) of normal-liquid '^He when plotted 
as a function of temperature for different molar volumes 
[5^. Rather sharp crossing points are also observed [9 in 
the specific heat curves c{T, X) of many heavy-fermion 
compounds, with X as a second thermodynamic variable, 
e.g., when measured at different pressures {X — P) as 
in CeAla [T^, or for different magnetic fields {X = H) 
as in CeCu6-a;Al:r HU, RuSr2Gdi.5Ceo.5Cu20io-5 [H], 
Mni„:rFe^Si [B], Mni_^Co^Si [TS], CuaOSeOg [S], and 



MnSi [ig. 

Crossing points of specific heat curves have been found 
in theoretical investigations of lattice models for cor- 
related electrons such as the one-band Hubbard model 
|9ll6l - fT5] . At half-filling the curves obtained by plotting 
c(T, U) vs. T for different values of the local Coulomb re- 
pulsion U always cross at two temperatures, irrespective 
of the type of the lattice, as seen in the case of nearest- 
neighbor hopping in d = 1 fTWT] . d = 2 and 
d = 00 [T7], as well as for long-range hopping in d = 1 
|26| . In particular, the crossing point at high temper- 
atures was found to occur at a nearly universal value 
c* /ks — 0.34 |18I27| . Subsequently crossing points were 
noticed also in the charge susceptibility of the half-filled 
extended Hubbard model in d = 1 f2S] and the thermal 
conductance of strongly correlated quantum dots |29| . 

The crossing points discussed above all occur when the 
temperature dependence of some physical quantity such 
as the specific heat is plotted for different values of a 
parameter, e.g., pressure, magnetic field, or interaction 
strength. But crossing points have also been observed 
in experimental data and theoretical results for the fre- 
quency dependence of various response functions, e.g., 
(i) in the optical conductivity a{oj,X) of electron doped 
Nd2_5Ce5Cu04, with X = 5 as the doping level [6 , and 
other high-Tc materials [30] as well as in the pyrochlore- 
type molybdate family R2M02O7 (R=Nd, Sm, Eu, Gd) 
[3r, where X now denotes the ionic radius r associated 
with the different elements R; (ii) in the Raman response 
function xi'^iT) of the high-T^ material HgBa2Cu04+5 
plotted as a function of the Raman shift u for differ- 
ent temperatures [32j and (iii) in the dielectric func- 
tion e(a;,r) of the colossal magnetoresistance material 
LaMnOs [33] and the reflectivity R{uj,T) of the colos- 
sal dielectric constant material CaCu3Ti40i2 [35 when 
plotted as a function of the photon energy uj for different 
values of T. 

Sharp crossing points in the Raman response x(a;,T) 
have been found in theoretical investigations based on the 
Hubbard model [7"35] and the Falicov-Kimball model [551 
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[37]. The spectral function A{uj, U) of the Hubbard model 
computed within dynamical mean-field theory (DMFT) 
also shows a sharp crossing point |38J. 

Whenever isosbestic points occur in a function f{x,p) 
two separate questions arise [9,: 

(i) Why do the curves cross at all, and 

(ii) why is the crossing feature confined to a narrow 
region despite significant changes in the parameter 
P [35]? 

As discussed by one of us ^ the existence of crossing 
points in specific heat curves can be explained by a sum 
rule for the change of the entropy S{T, X) with respect 
to X in the limit of T — >■ oo. Furthermore, their sharp- 
ness can be linked to the smallness of the susceptibilities 
x(")(r, X) = of strongly correlated systems in 

an expansion of the specific heat c(T, X) with respect to 
X, where £^{T,X) is the conjugate variable to X |40| . 

In this paper we concentrate on the second question. 
In Sec. [n] we discuss the local nature of isosbestic points 
and derive a general expression for the leading parame- 
ter dependence of a function f{x,p) exhibiting isosbestic 
points. This yields a generalization of the two-fluid model 
[1] and provides a scheme for the analysis of isosbestic 
points in existing data. In Sec. |III| we then apply our 
approach to the isosbestic point in the optical conduc- 
tivity of the FaHcov-Kimball model gT] within DMFT 
p2 -4r. Using the exact solution for this model the isos- 
bestic point is analyzed in detail regarding sharpness and 
existence. In Sec. |IV| we investigate isosbestic points in 
the spectral function A{uj,U) of the Hubbard model in 
DMFT. In Sec. |V] we apply our scheme to existing ex- 
perimental data to demonstrate its wide range of appli- 
cability. For this purpose we investigate the sharpness 
of isosbestic points in the Raman response x{'-^tT) of 
the cuprate compound HgBa2Cu04+5, in photoemission 
spectra / (w, T) of VO2 thin films, and in the reflectivity 
Riu},T) of CaCu3Ti40i2. 



II. PROPERTIES OF EXACT AND 
APPROXIMATE ISOSBESTIC POINTS 

A. Existence of crossing points 

A family of non-monotonic curves, obtained by plot- 
ting a quantity f(x,p) as a function of x for different 
values of a parameter p, will in general intersect. The 
crossing points are located along a curve x*{p) defined 



by M 



df{x,p) 



dp 



(1) 



If these points of intersection are confined to a narrow re- 
gion (in which case we refer to them as isosbestic points) , 
the value x*{p) depends only weakly on p. If the curves 



intersect at a single point, x* does not depend on p at 
all. 

Eq. ([T|) immediately explains why the absorbance of 
a mixture of two liquid solutions with individual con- 
centrations ci , C2 , whose total concentration is constant 
(ci +C2 — c), has a sharp isosbestic point |2]. Namely, in 
this case the absorbance a{uj,Ci) depends only linearly 
on the concentration ci, i.e., has the special form 



a(w,ci) = ciai(w) + (c - ci)a2(a;). 



(2) 



If ai and a2 coincide at some frequency uj* , i.e., ai{uj*) 
02(0;*), then 



da{u!, ci) 



5c 1 



0. 



(3) 



This implies that for all concentrations ci the absorbance 
curves intersect exactly at one frequency uj* , or the equiv- 
alent wavelength. Quite generally, whenever a system is 
a superposition of two components, where the sum of the 
densities is conserved, isosbestic points in the curves plot- 
ted for different densities are bound to occur |38I48| . In- 
deed, this argument is used by Uchida et al. [5] and Kezs- 
marki et al. ^T] to explain the well-defined isosbestic 
points in the optical conductivity 17(0;,^) of numerous 
correlated electron materials, where p is some control pa- 
rameter, e.g., the temperature or the bandwidth. Start- 
ing from the sum rule for the optical conductivity, and 
assuming (i) that the spectral weight is only redistributed 
between electrons with two different energy scales, and 
that (ii) the optical spectrum can be decomposed and lin- 
early interpolated by two terms as a function of p, they 
arrive at an expression for the optical conductivity given 
by (t(cj,p) = ptTi(w) -|- (1 — p)(T2(w), which is identical to 
the result for the absorbance in Eq. ^ and hence leads 
to an isosbestic point at some particular frequency tj* 

m- 

Before we continue with our discussion we need to in- 
troduce a suitable nomenclature. Based on the behavior 
of the curve of crossing points x*{p) we identify three dif- 
ferent cases: (i) A globally exact isosbestic point, which 
is characterized by x*{p) — const, corresponding to a 
complete p-independence of f{x*,p): (ii) a locally exact 
isosbestic point, which arises when x*{p) exhibits a local 
extremum or a higher order stationary point around some 
value p — po, corresponding to a locally p-independent 
f{x*{p),p) around p — Po', (hi) an approximate isos- 
bestic point, which corresponds to a weak p-dependence 
of x*{p). 

We note that locally exact isosbestic points can be fur- 
ther classified by the order of the stationary point of x* (p) 
at pq. If 



d'^x*{p) 



dp*' 



= 



(4) 



for k < n, the isosbestic point is called an "isosbestic 
point of nth order". 
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B. Sharpness of isosbestic points 

To understand the general origin of sharp isosbestic 
points in a function f{x,p), it is important to note that 
isosbestic behavior is usually observed only in a certain 
parameter range around some particular value po and 
can be expected to break down away from pq. Accord- 
ingly, the weak p-dependence of x*{p) required for an 
approximate isosbestic point will generally be a local 
phenomenon, i.e., it applies only for a finite p-interval 
around pq. This observation of the local nature of isos- 
bestic points motivates the following expansion around 
Pq 

f{x,p) = f{x,po) + {p- Pa) Fi{x,po) + O [(jp-pof] , 

(5) 

where Fn{x,p) = One has Fi{x*{po),po) 

= 0, by virtue of the definition of x*{p) [Eq. (fTl)]. As a 
consequence, the linear approximation of f{x,p) 

/approx(a;, p) = f{x,Po) + {p- Po) Fi {x, po) , (6) 

will exhibit an exact isosbestic point at a; = x*{pa) for 
all values of p. Since 



f{x,p) = /approx(a;,p) +0 [{p-po)'^] 



(7) 



the second term is seen to be responsible for deviations 
from the exact isosbestic point described by /approx(a^,p)- 
This result can be put to a test by checking the validity 
of the following approximation: 

f{x,p) = f{x,p)-{p-po)F,ix,po) (8) 
= f{x,p„)+0[ip-po)^] (9) 
^ f{x,Po). (10) 

The quantity f(x,p) represents the reference curve 
f{x,po) plus all terms which are responsible for the de- 
viations from the exact isosbestic case. The verification 
of its approximate parameter independence, i.e., the ver- 
ification of the approximation f{x,p) ~ f{x,po) allows 
one to quantify the importance of the subleading terms 
O [{p — Pn)^] and therefore the deviations from the exact 
isosbestic case. 

In general, the value of po should correspond to the po- 
sition of the weakest p-dependence of x*{p) and, without 
any further information on the system, is generally hard 
to determine unambiguously. Of particular interest how- 
ever are the cases where po corresponds to an extremum 
or some higher-order stationary point of x*{p). These 
cases lead to locally exact isosbestic points where po is 
determined by x*' {po) — 0. The new variable {p — po) 
is then a small parameter of the system, e.g., an inter- 
nal frequency scale if p is a frequency. From the implicit 
Eq. one obtains 



dp 



c*{p) = - 



dFi{x,p)/dp 
dFi{x,p)/dx 



(11) 
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Figure 1: Isosbestic point in the optical conductivity of 
the Falicov-Kimball model at T = l4l] (Gaussian DOS 
[Uc — 1.5] with density n of mobile and density wi = 1/2 
of immobile particles). Spectral weight is transferred from 
low frequencies to the Hubbard peak as electron-density is in- 
creased, producing an isosbestic point at medium frequencies. 



which means that for x*'{po) — the subleading terms in 
Eq. [rjare at least of order (p — po)^, thus corresponding 
to particularly weak deviations from the exact isosbestic 
point in /approx(2;,p) in a vicinity of po- In practice isos- 
bestic points corresponding to an extremum of x*{p) are 
rather common, as will be discussed in the following sec- 
tions. A more mathematical treatment regarding the be- 
havior around such stationary points of x* (p) is provided 
in Appendix [K\ 

In the rest of the paper we will use these results to show 
that isosbestic points are indeed connected to a leading 
parameter dependence as described by Eq. (6|. 

Finally we emphasize that it is not possible to quan- 
tify the sharpness of isosbestic points beyond Eqs. (8p0 I, 
since the perceived sharpness always depends on the mag- 
nification of the crossing region, which is not an intrinsic 
property of f{x,p) itself (see Appendix IbI) . 



III. ISOSBESTIC POINTS IN THE OPTICAL 
CONDUCTIVITY OF THE FALICOV-KIMBALL 
MODEL 

A. Model and optical conductivity in DMFT 



The Falicov-Kimball model |41I50| describes itinerant 
electrons (c-electrons) interacting with localized electrons 
(/-electrons) by a local Coulomb repulsion. The Hamil- 
tonian is given by 



H 



(12) 



// fi 



x—x*{p) 



It can be viewed as a Hubbard model with spin- 
dependent hopping and two different chemical potentials. 
The symbols cj (fj) and Ci (fi) denote the itinerant 
(localized) electron creation and annihilation operators. 
The chemical potential fi constrains the total number of 
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c-electrons, while Ef is an orbital energy of the localized 
electrons. Here U is an on-site repulsion between the two 
electron species. 

At low temperatures, the optical conductivity of the 
Falicov-Kiniball model with the Gaussian density of 
states (DOS) is known to exhibit a distinct isosbestic 
point [41J when plotted for various c-electron densities 
(Fig. [T|. In the following we will discuss the existence 
and sharpness of isosbestic points and concentrate on the 
ungapped system at zero temperature and use the semi- 
elliptic DOS p{e) = (27rt)- V4t2 - with t = 1. 

Within dynamical mean-field theory (DMFT), the 
model can be solved exactly for the Green function G{uj) 
or the self-energy 5](a;) owing to the simplicity of the 
impurity problem |51| : 



G(w) 



I — Wi 



Wi 



(13) 



where GqI{co) = S(ci;) + G~^{lu) denotes the Green func- 
tion of the effective medium (Weiss field). The weight wi 
equals the average /-electron concentration which will be 
set to 1/2 in the following. Due to particle-hole symme- 
try, it is sufficient to restrict the treatment to c-electron 
densities n < 1/2, i.e., < 0. 

The optical conductivity in DMFT is essentially given 
by the particle-hole bubble [52^, which on the real axis 
corresponds to 

/oo r 
dep{e)v{ef / duj' (w') {J + oj) 
-OO J 

[J) ~ J {lo' + io) 



A^{u)) = — I/ttIui [cj + — I](iuj) — e] 



(14) 
(15) 



where A^[lo) denotes the spectral density and p{e) the 
density of states (DOS). We fix the velocity function u(e) 
according to the treatment in Ref. j 53| . 



B. Existence of isosbestic points in (7{u)) 

The formation of isosbestic points in the optical con- 
ductivity (T(a;) when plotted for different doping levels of 
the c-electrons in Fig.[T]can easily be deduced on physical 
grounds, viewing the Falicov-Kimball model as a disor- 
dered medium in which spinless electrons interact with 
localized impurities. In this picture the scattering rate 
T~^(/i) oc — ImSp(aj = 0) is monotonically related to the 
c-electron density through the chemical potential /i, and 
reaches a maximum at half-filling (see Fig. [2|d). Since 
r~^(0) diverges at the metal-to-insulator transition, it is 
always possible to suppress the DC conductivity of the 
half-filled FK model simply by increasing U . Thus, in the 
intermediate ?7-regime aociU, n) will have a minimum or 
a dip around half-filling which means danc/dn < 

Besides the mathematical rea- 
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Figure 2: Doping dependence of (a) the rescaled DC con- 
ductivities U^aoc (p) and (b) self-energies ImE^(0)/(7^ for 
the semi-elliptic DOS. With growing interaction ImE^(O) be- 
comes increasingly peaked with an extremum at half-filling 
{/i = 0). Therefore the scattering rate oc — ImE(O) be- 
comes maximal at half-filling. The metal-insulator transition 
at [7 = 2 is distinguished by a vanishing DC conductivity 
and a diverging scattering rate t~^. For sufficiently high, 
yet finite U this generic property of the FK model leads to 
the formation of minima in aoc (a*) around for jj. — 0. This 
corresponds to dcTDcifJ-) /d^ < 0. Note that there exists a 
one-to-one mapping between p and the carrier density n. In 
Appendix [e| we show that S;j(a;) — 'S{uj + p). 



for non-zero frequencies and for sufficiently strong U 
can be established by taking into account the formation 
of the well-known charge-transfer or Hubbard peak at 
LUc{U,n) ~ U (see Fig. [TJ, which is a property of the 
optical conductivity of correlated systems. Representing 
a robust indicator of the increasing correlation between 
/- and c-electrons, the peak becomes increasingly pro- 
nounced as the system is doped towards half-filling. For 
sufficiently strong U we thus have da (wc) /dn > and 
dcTDc/dn < 0, and as a consequence there exists a cross- 
ing point uj*{n,U) (c.f. Eq. Q) with < uj* < uc and 
da (w*, n) /dn = as depicted in Fig. [l] The formation 
of isosbestic points is thus an effect of the correlation as- 
sociated transfer of spectral weight from low to higher 
frequencies as the system is doped towards half-filling. 
The argument does not depend on the specifics of the 
density of states under consideration and thus applies to 
Fig. [T| (Gaussian DOS) as well as to Fig. [s] (semi-elliptic 
DOSV 



C. Sharpness of isosbestic points in (j{uj) 

The symmetric density dependence of the crossing fre- 
quency (jj*{n) depicted in Fig. |4] reflects the important 
role of particle-hole symmetry for locally exact isosbestic 
points in a {uj,n). Because of the existence of the one-to- 
one correspondence between the chemical potential p and 
the carrier density n in the ungapped system both param- 
eters are equivalent. The chemical potential, however. 
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Figure 3: Optical conductivity (a) a{uj, fi (n)) and (b) the 
quantity a{ij, jj, (n)) for the Falicov-Kimball model with the 
semi-elliptic DOS. The curves cr(tj,/i(n)) in (b) essentially 
lie on top of a (iu, jj, — 0) for the substantial density inter- 
val n G [0.35 ... 0.5]. Thus, the parameter dependence in 
a{Lj, n (n)) essentially stems from the first contribution, which 
is accompanied by an isosbestic point. The second derivative 
of (T (oj, /i) w.r.t. /i, which is needed for the calculation of 
a{u},ii), is obtained from the exact solution |54| . 



represents the more convenient parameter because of the 
simple /z-dependence of the self-energy S(a;) (see Ap- 
pendix[E| and of the spectral densities A^(uj) in Eq. (14 1. 
Thus, concentrating on ungapped systems, we expand 
around half-filling, i.e. /x = 0, taking fx as a. small param- 
eter 



a (lj, y)^a {u, 0) + {u, ^l = 0)+O {^i*) . (16) 

With Eq. ^ we are directly led to the conclusion that the 
first two terms constitute an exact isosbestic point, while 
corrections of order O (/z*) introduce the /i dependence 
into uj* (/i) (or a;*(ri)). According to Eq. ([s]) we therefore 
introduce 

1 52 

= a{Lu,^i^o) + o[^i% 

and test its approximate /i independence in a vicinity of 
oj* (/i = 0). This is confirmed graphically in Fig. |3] the 
obvious agreement with the reference curve can be seen 
as a measure of the weak influence of higher order con- 
tributions. This provides an especially simple and mean- 
ingful interpretation of the approximate isosbestic points 
in a{iJ,n): Namely, plots of the conductivity showing 
this behavior exhibit a density dependence which can be 
treated in perturbation theory for the parameters under 
consideration. 
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n 

Figure 4: Isosbestic points uj* (n, U) of the ungapped sys- 
tem using the semi-elliptic density of states. The symmetry 
around half-filling is a consequence of particle-hole symmetry. 



IV. ISOSBESTIC POINTS IN THE SPECTRAL 
FUNCTION OF THE HUBBARD MODEL 

We now apply our approach to the prominent isos- 
bestic point in the spectral function A (w, U) of the 
particle-hole symmetric Hubbard model within DMFT 
for which, employing symmetry arguments, the crossing 
frequency was shown |38| to obey the following depen- 
dence on U, 



J* {U) = Co + ciU^ + 0[U^ 



(17) 



For a symmetric DOS with van-Hove singularities at the 
band-edges (at ±5) like the semi-elliptic DOS, Eq. (17 1 
takes the more specific form oj* (U) — b + ciU^ + 0{U^. 

From the infinite slope of the spectral function A{uj,0) 
at w = cj* (0) it then follows that, in the limit U 0, 
one can never obtain a sharp crossing point defined by 
the coordinates {x{U),y{U)) with x{U) = lo* {U) and 
y{U) — A{uj*{U),U). In order to avoid this addi- 
tional mathematical complication, let us here consider 
the Gaussian DOS with infinite bandwidth, since this 
DOS is a smooth function so that a weak parameter de- 
pendence ofuj*{U) always translates to a sharp isosbestic 
point. From particle-hole symmetry, we have the follow- 
ing perturbative expansion 



A{oj,U) 



n>0 



C/'"^2„(C^) 



with 



A2iuj) = --^Imly 



dep{e) 



E2(w) 



[uj + iO+ -eY 



(18) 



(19) 



Here J^2{^) denotes the self-energy in second order per- 
turbation theory [55 in U. We thus obtain from Eqs. ^ 

lol 



A{uj,U) = A{lu,U) -U^A2iuj) 



(20) 



and again obtain A {uj, [/) ~ A {u, 0), as can be seen from 
Fig-i 
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Figure 5: Numerical renorinalization group calculation [56157] 
of the spectral function for the half-filled Hubbard model 
(a) using the Gaussian DOS p(e) = 4/^exp[-'(4e)^]. The 
quantity A {uj, U) shows practically no fZ-dependency up to 
Ui ~ 0.35 (b). The deviations from the locally exact isos- 
bestic point at higher U stem from higher order perturbations 
and become noticeable only above Ui. 



V. ANALYSIS OF EXPERIMENTALLY 
OBSERVED ISOSBESTIC POINTS 

The generality of Eq. ([7| allows us to investigate the 
sharpness of a whole class of isosbestic points in entirely 
unrelated systems within a common framework. An ex- 
isting set of theoretical or experimental data of a quantity 
Q{x,pi) obtained for a number of different parameters, 
i = 1, 2, ... n can be analyzed in the vicinity of an isos- 
bestic point x*{p) by comparing it to 

Qapprox = Q{x,Po) + {P~ Po)Ql{x,Po)- (21) 

Here the function Qi {x,pq) is approximated as 
Q{x,Pa) - Q{x,pi3) 



i{x,pq) 



{Pa - Po) - {PfS - Po) ' 



(22) 



with PonPp G [pi---Pn] sufficiently close to po to ob- 
tain a proper approximation of Qi{x,pq). As before the 
applicability of Eq. (21 1 is then tested by verifying the 



approximate p independence of 

Q{x,pi) = Q{x,pi) - (p- 
- Q{x,Pq). 



Pq)Qi{x,Po) 



(23) 
(24) 
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Figure 6: (a) Antinodal Raman response x ("^i T) of 
HgBa2Cu04+5 [32] at optimal doping and (b) the function 
x{^,T). The weak temperature dependence of x(aj,T) ex- 
plains the sharpness of the isosbestic point and also confirms 
the validity of ansatz ( 25 1 within the available data. 



HgBa2Cu04+5 for various temperatures at optimal dop- 
ing ((5 = 0.16); which shows a distinct isosbestic point 
for up to ~ 100/^. Proceeding as described above we 
make the following Sommerfeld-type ansatz for the low 
temperature dependence of x {^i T): 

X {u, T) = X (w, 0) + T\2{^) + 0[T\ 



where 



X2(t^) 



^1 



(25) 



(26) 



as in Eq. (22). Again, the validity of Eq. (25) can be 
tested via the T-independence of the quantity 

xi^-X) = x{^.T)-T'^xiiy:) 

~ X(c^,r = 0). (27) 

This is done in Fig. |6] which shows good agreement with 
Eq. (27). To determine X2{^) 

we chose T\ — \AK and 
T2 = 90iir, which is by no means a unique combination 



It is interesting to note that the isosbestic point 
thus the quadratic T dependence in Eq 



(25) 



and 
are specific 

features of the optimally doped system and are absent 
for other dopings [35] (at least for the provided tempera- 
tures). For the case of optimal doping, we have explained 
the sharp isosbestic point in x(a;,r) by its essentially 
quadratic temperature dependence. 



A. Raman response of HgBa2Cu04+5 

We now apply the scheme to the isosbestic point 
observed in measurements |32| of the antinodal (-Big) 
Raman response %(a;,T) of the cuprate compound 



B. Photoemission spectra of V02/Ti02 (001) thin 
films 

Next we apply the method to explain photoemission 
spectra /(w,T) of V02/Ti02 thin films [59 , which ex- 
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Figure 7: Temperature dependence of (a) photoemission spec- 
tra I {u},T) of a VO2 thin film [59| . and (b) the function 
I {ui,T) defined in Eq. (29 1. The weak temperature depen- 



dence of J (cj, T) explains the sharpness of the isosbestic point 
and also supports the importance of electron-phonon interac- 
tion. 
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Figure 8: Temperature dependence of (a) the reflectivity 
R{ijJ,T) of CaCu3Ti40i2 [Ml and (b) the function R{uj,T). 
The weak temperature dependence of R {uj, T) explains the 
origin of the isosbestic point. 



hibit a marked isosbestic (see Fig. [7|. The strong tem- 
perature dependence of / (w, T) was seen |59) as an indi- 
cation of strong electron-phonon coupling. Consequently, 
I {uj,T) was jM] convincingly reproduced through the 
spectral function A (w, T) of the independent boson 
model through I {uj,T) cx A(uj,T). The temperature de- 
pendence of A (w, T) with A(a;,r) = A{uj,0)+TAi {uj) + 
O (T^) justifies the following ansatz: 



/(w,r) = /(w,0) +T/i {ij)+0[T^ 



(28) 



where 



Ji(^) 



IiL0,T,)-I{LU,T2) 

TI-T2 



is extracted from the data (Ti = 150K and T2 — 
200K) . The ansatz ( 28 1 is confirmed by the relative T- 
independence (well within measurement accuracy) of the 
quantity 



i{lo,t) = /(c^,r)-r/i(w) 
~ /(c^,r = o), 



(29) 



plotted in Fig. [TJd. This provides an explanation for the 
sharpness of the isosbestic point, and the leading lin- 
ear temperature dependence supports the importance of 
bosonic excitations along the line of Ref . [59 . 



C. Phonon modes in CaCu3Ti40i2 

We investigate isosbestic points in refiectivity R{uj,T) 
measurements of CaCu3Ti40i2 around phononic excita- 
tions 534] . Again, assuming the leading T-dependence 



to be linear, we proceed along the lines of the preceding 
section and introduce 



i?(w,T) =i?(w,T)-Ti?i( 



where 



Ri (w) 



R{io,Ti) - R{lu,T2) 
T1-T2 



(30) 



is extracted directly from data. Note that the sensitiv- 



ity of Eq. ( 30 1 on the quality of the available data re- 



quired the use of a comparably large temperature dif- 
ference (Ti = 4QK, T2 = 175 A'). The high quality of 
the match in Fig. |8] explains the sharpness of the isos- 
bestic point and also shows that a linear temperature 
dependence already provides a proper description of the 
available measurement data. 



VI. CONCLUSION 

In this paper we developed a framework which allows 
one to extract information about a correlated system 
by analyzing the sharpness of an approximate isosbestic 
point in a physical quantity f(x,p). Our central result 
is given by Eq. ([t]) which corresponds to a straightfor- 
ward generalization of Eq. Q , describing exact isosbestic 
points [1'2]. We established a direct connection between 
the sharpness of isosbestic points and the existence of a 
small parameter Sp ^ p — po, which can be used to de- 
scribe the leading parameter dependence of the system in 
the parameter interval under consideration. In the case 
of additional symmetries pq can provide useful informa- 
tion about the system, e.g., about relevant energy scales. 
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Furthermore, symmetries can lead to a high power m ap- 
pearing in the correction term leading to locally 
exact isosbestic points as found in some of the examples. 

Experimental or theoretical insight regarding pq, m, or 
even the shape of Qi{x,po) [Eq. p2| )] is straightforward 
to test against measured data. We perform this kind of 
analysis for the Raman response of HgBa2Cu04+5. For 



this system, we suggest a reinvestigation of (25l-(29l us 



ing higher resolved data (especially for lower tempera- 
tures), in order to determine the function X2(w) better 
[Eq. (25)]. In particular it has a simple temperature de- 



pendence only for optimal doping, which raises the ques- 
tion of whether other dynamical quantities exhibit it as 
well. Another interesting question is whether other (opti- 
mally doped) cuprate compounds have similar isosbestic 
points in their Raman spectra. 
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Appendix A: Size of the crossing region around 
locally exact isosbestic points 

We quantify the influence of the higher-order contribu- 
tions, i.e., the terms responsible for the deviations from 
x*{p) = const. For this we investigate the p-dependence 
of an isosbestic point at {x*{p),p), with x*{p) given im- 
plicitly by Eq. ([ij . We will estimate the size of the region 
where the curves of f{x,p) cross, e.g., for 



X = x*{po) + Sx . 
P = Po + Sp, 



(Ala) 
(Alb) 

in the vicinity of a certain parameter value pq. 

We first approximate f(x,p) near po by a function that 
is linear in p, 



/approx(a;,p) f{x,Po) + SpFi{x,po) 



(A2) 



where Fn{x,p) — cf^ f{x,p) / dp'^ . Equation (A2| is exact 
at Pq: /approx(a;,po) = /(a;,Po)- As in the two-fluid model 
that was discussed in the introduction, /approx(a;,p) has 
an exact isosbestic point at (a;*(po)jPo) because i^i(a;,po) 
vanishes there. 

Next we determine the difference between the function 
f{x,p) and the linear approximation /approx(2^7P) is- For 
this purpose we flrst expand Fi{x,p) near x*{p) for some 
fixed p and obtain 



Fi{x,p) = Fiix*ip),p) 



^ [x~x*{p)]" d'^Fi 



n=l 



dx^' 



x—x* (p) 

(A3) 



By virtue of the definition of x*{p) [Eq. ([T|)[, the x- 
independent term vanishes, i.e., Fi{x*{p)^p) — 0, leading 
to 

i^i(x,p) = c[p) [x - x-'ip)] + 0{[x - x*{p)f) , (A4) 

with c{p) = dFi{x,p)ldx\^=^*(py 

The height of the crossing region is thus given by 

Sf{x,p) = f{x,p) - /approx(a;,p) 



Sp^ 



F2{x,po) + ^Fs{x,po) + 0{p'). (A5) 





Using Eq. ( A4 1 we can now express the higher derivatives 

as (n > 0) 

i^„+i(x,p) = c(")(p) [x-x*{p)] 

'E(I)^^""'^(^)-'^''H^^)' (A6) 

where c(")(p) = d"Fi{x,p)/d^x\,=,.(^py 

So far the parameter value p has not been specified. 
We now suppose that x*{p) has an extremum at po, in 
which case we term (x*(po)iPo) a first-order isosbestic 
point. From x*'{po) = it then follows that 

Sf = Sp"" Sx + OiSp^Sx"^) . (A7) 

If instead x*{p) has a saddle point at po (second-order 
isosbestic point), then furthermore x* (po) — 0, yielding 



Sf 



c'{p) 



Sp^ 



clip) 
6 



V Sx + 0((5p*(5a;") , (A8) 



and similar for higher-order isosbestic points. We thus 
conclude that an nth-order isosbestic point, for which 
the first n derivatives of x*{p) vanish at po, will have a 



narrow crossing region 



ri+l 



Sf = 5xY^ c^^-^\p) 5p^ + OiSp'^+^Sx") , (A9) 



fc=2 
„n+2 



0((5p"+") if 5x = 



(AlO) 



i.e., the first n+1 powers of Sp are suppressed by a factor 
Sx, and the contribution of order Sx^ is suppressed by a 
factor Sp"^'^. The last equation applies if we consider 
the height of the crossing region Sf{x*{po),p) for various 
parameters p near po, and this height is only of order 
0(5p"+^) for an nth order crossing point. We conclude 
that the isosbestic point x*{p) of a function f{x,p) will 
be particularly sharp in the vicinity of a extremum or 
higher-order stationary point of x*{p). 



Appendix B: Quantification of the sharpness of 
isosbestic points 

It is not possible to introduce an unambiguous mea- 
sure for the sharpness of an approximate isosbestic point 
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because of its explicit dependence on the magnification of 
the crossing region. This means that the observed sharp- 
ness always depends on the specifics of plotting. For fixed 
plot parameters, i.e., for f{x,p) plotted for x € [xa,Xb\ 
and p g [pi,pn], one can, however, define a measure for 
the relative sharpness of a crossing region through the 
ratio between the area of the plot and the area of the 
crossing region. For this we define the vertical width of 
the plot as 



A/ = max \f{xi,pi 



f{x2,Pi)\ 



(Bl) 



with i = 1, . . . ,n and Xi^2 G [xa, xt] and correspondingly 
the vertical width of the crossing region as 



A/* =max /* 



(B2) 



where /* = f {x* (pi) , pi) . With the horizontal width of 
the crossing region defined by 



Ax* = max \x*{pt) - x*{pj)\ 
id 

one can introduce a quantity S through 



S = 



Af* Ax* 



(B3) 



(B4) 



where Ax = \xa — Xi,\ which quantifies the sharpness, 
i.e., S — oo for an exact isosbestic points and S ^ I 
for sufficiently sharp isosbestic points. We stress the ex- 
plicit dependence of S on Ax and the range of parame- 
ters [pi,p„] unless S = oo. For finite S the qualification 
"sharp" should thus be considered a relative term and be 
used with caution. 



Appendix C: Existence of isosbestic points in u{uj) 

First we investigate the system specific mechanisms 
that lead to isosbestic points in cr (w, n). To allow for a 
rigorous analysis, we restrict ourselves only to cases in 
which there exists at most one isosbestic point. This ap- 
plies to the results [H] obtained for the Gaussian DOS as 
well as to our calculations for the semi-elliptic DOS. Due 
to particle-hole symmetry one has a {oJ,n) = cr {uj, —fJ.). 
According to Eq. Q, isosbestic points are determined by 



d(j{uj*{n), n) 



dn 



(CI) 



u}* (n) 



In view of the one-to-one mapping between the chemical 
potential fi and the carrier density n in the ungapped 
system, which is a direct consequence of the Green func- 
tion Giuj, n) being a function only of Lu + fi (see Appendix 
Ie|), this is equivalent to 



0. 



(C2) 
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Figure 9: DOS with finite support {jj. < 0, the shaded area in- 
dicates filled states): The term involving (— w) in Eq. (C4| 
vanishes for all frequencies oj > wq, while (oi) > for all 
(jjo < oj < ioi . Thus, the derivative da {luq) /dn will necessary 
be positive for all cjo < w < lji; for even higher frequencies it 
vanishes identically. 



Since S (w) is also only a function of io + ^ (see Appendix 
IeI), we obtain 



^g(^, m) 
dfi 



d 



= ao I dep{e)v{ey / du' — [A, {oj') A, {oj' + 



f{Lu')-fiuj'+i0) 



(C3) 



Integration by parts and using the property df{io)/dLU = 
—6(uj) of the Fermi function gives 



d ao 
—a{uj,fi) = — 



dep{e)A, (0) [A,{u)-A, (-c^)] . 

(C4) 

at T = 0. 

In Appendix [D] we show that for less than half- filling, 
the derivative da{uj, p) / dp, is positive around a > 
and never becomes negative in the high frequency limit. 
Since we assume the existence of at most one isosbestic 
point, it follows that it is exclusively determined by 
the /i-dependence of the DC conductivity udc (m) = 
a{uj^ 0,/i): 



d_ 
dp 



(C5) 



for less than half- filling. This relation expresses the 
equivalence between the existence of isosbestic points 
and the formation of a suppression of DC conductivity 
a DC ip) resp. ajjc {n) around half-filling. 



Appendix D: high-frequency limit 

We determine the sign of d(T{ujo)/dp in the high- 
frequency limit, taking p < 0. In this case, the Green 
function and the self-energy are shifted to higher ener- 
gies by p compared to the half- filled case. For a symmet- 
ric DOS with finite bandwidth it follows from Eq. (C4| 



that there exists a frequency ujq > with A{ljo) > and 
A{-ujo) = (c.f. Fig. |3|. This yields da{uJo)/dp > 
because the contribution involving A^(—ujq) in Eq. (C4| 
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must identically vanish due to the finite support of p(e) 
and p(e). Symmetric densities of states with infinite 
bandwidth behave similarly: There exists a ujq with 
A (uj) > A (— w) for all uj > uq because of the stronger 
suppression of the contribution involving A (— cj). This 
analogously leads to da{uJo)/dfJ. > 0. 



Appendix E: Shape invariance of single-particle 
quantities with respect to doping 



representing the first few terms of its skeleton expansion. 



As a consequence of Eq. ( 13 1 , the model has an especially 



simple doping dependence, where a change of fj, (for con- 
stant wi) merely shifts the Green function and the self- 
energy, leaving their shape invariant. Mathematically, 
this is expressed by the fact that oj and fi enter G{uj) and 
I](a;) only through uj + ^. To see that, we consider the 
inverse [46 of the Hilbert transform R [G(a;)], expressing 
Lu + fi ~ 'E{uj) exclusively through G{uj). Equation (13 1 
then takes the form 



Within DMFT, the self-energy of the model can be ex- 
pressed explicitly in terms of the interacting Green func- 
tion |41|51| (wo = 1 - wi): 



U 
2" 




G(a;) 



AwqU 
GM 



(El) 

representing the exact summation of the skeleton expan- 
sion for S(aj). In the low-U regime an expansion of S](w) 
in terms of the small parameter UG{uj) gives {wi = 1/2) 

S(c.) ^lu+l U'Giu;) - 1 U^G^u;) + 0[U% (E2) 



Giu) 



1 — Wi 



UJ + H- R[G{uj)]+G-^uj) 

wi 

u + fi-R [G{uj)] + G-i(w) - U ' 



(E3) 



implicitly defining the Green function as function of w-h/z, 
i.e., G{uj) = F{ijj + For the special case of a semi- 



elliptic DOS, Eq. (E3l leads to a cubic equation, which 
can be solved explicitly [54J for G {uj + fi). 
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